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We have reexamined the possibihty of explaining the solar neutrino problem through long- 
wavelength neutrino oscillations induced by a tiny breakdown of the weak equivalence principle 
of general relativity. We found that such gravitationally induced oscillations can provide a viable 
solution to the solar neutrino problem. 



Nature seems to be most strongly in agreement with 
neutrino oscillations. The compelhng evidences coming 
from solar neutrino experiments , that span over two 
decades, and from atmospheric neutrino experiments |^] 
are difficult, if not impossible, to be accommodated with- 
out admitting neutrino flavor conversion. Nevertheless 
the dynamics underlying such conversion is yet to be es- 
tabhshed and in particular does not have to be a priori 
related to the electroweak force. 

The interesting idea that gravitational forces may in- 
duce neutrino mixing and flavor oscillations if the weak 
equivalence principle of general relativity is violated, was 
proposed by Gasperini |^ and independently by Halprin 
and Leung Q about a decade ago, and thereafter, many 
works have been performed on this subject ||^-[l^. In 
Ref. |l5| this was shown to be phenomenologically equiv- 
alent to velocity oscillations of neutrinos due to a possi- 
ble violation of Lorcntz invariance [^6|. So even a tiny 
breakdown of the space-time structure of special and/or 
general relativity may lead to flavor oscillations even if 
neutrinos are strictly massless. 

Some theoretical insight on the type of gravitational 
potential that could violate the weak equivalence princi- 
ple can be found in Ref. |jl^ . A discussion on the depar- 
ture from exact Lorentz invariance in the standard model 
Lagrangian in a perturbative framework is developed in 
Ref. @. 

Several authors have investigated the possibility of 
solving the solar neutrino problem (SNP) by such gravi- 
tationally induced neutrino oscillations P-pl|, generally 
flnding it necessary, in this context, to invoke the MSW 
like resonance since they conclude that it is impos- 
sible that this type of long-wavelength vacuum oscilla- 
tion could explain the specific energy dependence of the 
data |jl^. 

Recently these neutrino oscillation mechanisms have 
been investigated ||l^-[l^ in the light of the experimental 
results from Super-Kamiokande (SK) on the atmospheric 
neutrino anomaly, obtaining stringent limits for the i^^ 
v-r channel. 

We consider in this letter the possibility of explaining 



the most precise and recent solar neutrino data coming 
from gallium, chlorine and water Cherenkov detectors by 
means of neutrino mixing due to a "just-so" violation of 
the weak equivalence principle (VEP). We demonstrate 
that all the data can be well accounted for by the VEP 
induced long-wavelength neutrino oscillation in contrast 
to previous conclusions P, p^|jT5| ]. 

We assume that neutrinos of different species will in- 
cur different time delay due to the weak, static gravita- 
tional field in the intervening space on their way from 
the Sun to the Earth. Their motion in this gravitational 
field can be appropriately described by the parametrized 
post-Newtonian formalism |]l9| with a different parame- 
ter for each neutrino type. In this manner neutrinos that 
are weak interaction eigenstates and neutrinos that are 
gravity eigenstates will be related by a unitary transfor- 
mation that can be parameterized, assuming only two 
neutrino flavors, by a single parameter, the mixing angle 
6*13 which can lead to flavour oscillation . 

Let us briefly revise the formalism that will be used 
in this work. We will assume oscillations only between 
two species of neutrinos, which are degenerate in mass, 
either between active and active {u^ <-> v^,Vt) or active 
and sterile {v^, <-> Vg, Ug being an electroweak singlet) 
neutrinos. 

The evolution equation for neutrino flavors a and (3 
propagating through the gravitational potential 4>{r) in 
the absence of matter is [|j : 



dt 



E(t>{r)A-f 



cos 29q 


sin 29q 






sin 29g 


— cos 29g 







(1) 



where E is the neutrino energy; A7 is the quantity which 
measures the magnitude of VEP, it is the difference of 
the gravitational couplings between the two neutrinos in- 
volved normalized by the sum. 

There are many possible sources for cf), but it is gen- 
erally believed that the Super Cluster contribution (0 ~ 
3 X 10~^) would be the dominant one |^0|. Therefore, 
it seems reasonable to ignore any variation of over the 
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whole solar system and take it as a constant . In this 
case Eq. (|l|) can be analytically solved to give the survival 
probability of v^, produced in the Sun after traveling the 
distance L to the Earth: 



1 -sin^ 26*0 sin^ — , 
A 



where the oscillation wavelength A is given by, 

MeV" 



A 



TT km 



5.07 



E 



(2) 



(3) 



which in contrast to the wavelength for mass induced 
neutrino oscillations in vacuum, is inversely proportional 
to the neutrino energy. 

In this case the survival probability is a function of two 
unknowns parameters that can be fitted, or constrained, 
by experimental data: A7 and sin 20(3. Since the value 
of the potential (j) in our solar system is somewhat un- 
certain we will adopt the procedure used by other 
authors and work with the product 0A7. 

We will perform a fit of the rates and SK recoil-electron 
spectrum but not take into account the day night effect 
(or zenith angle dependence) in SK. This is justified by 
the fact that day night variations can not be induced by 
this mechanism, and therefore, are irrelevant in deter- 
mining the allowed parameter region. We will comment 
about the possible seasonal variations at the end. 

We first examine the observed solar neutrino rates in 
the VEP framework. In order to do this we have cal- 
culated the theoretical predictions for gallium, chlorine 
and Super-Kamiokande water Cherenkov solar neutrino 
experiments, as a function of the two VEP parameters, 
using the solar neutrino fluxes predicted by the Stan- 
dard Solar Model by Bahcall and Pinsonneault (BP98 
SSM) taking into account the eccentricity of the 
Earth orbit around the Sun. 

We then have performed a analysis to fit these pa- 
rameters and an extra normalization factor fs for the 
®B neutrino flux, to the most recent experimental re- 
sults coming from Homestake ^ Rqi = 2.56 ±0.21 SNU, 
GALLEX g and SAGE || combined i?Ga = 72.5 ± 5.5 
SNU and SK g Rqk = 0.475±0.015 normalized to BP98 
SSM. The deflnition of the function to be minimized 
is the same as the one used in Ref. which essen- 
tially follows the prescription given in Ref. except 
that our theoretical estimatives were computed by con- 
voluting the survival probability given in Eq. with 
the absorption cross sections taken from Ref. [25[ and 
the neutrino-electron elastic scattering cross section with 
radiative corrections |^ and the solar neutrino flux cor- 
responding to each reaction, pp, pep, ^Be, ^B, ^'^N and 
^^O and other minor neutrino sources such as ^''F or hep 
neutrinos are neglected. 

We will first discuss our results for active to active 
conversion. We present in Fig. 1 (a) the allowed region 
determined only by the rates with free fg and in Table || 
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FIG. 1. Allowed region for sin^ 2dG and |(?!)A7| for (a) 
the rates only, (b) SK spectrum only and (c) rates -I- SK 
spectrum combined. The best fit points are indicated by the 
crosses and the local best fit points in the other 90 % C.L. 
islands are indicated, in each plot, by the open circles. The 
"test point" which will be used in Fig. 2 and 4 is indicated 
by the open square (see also the text). 
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FIG. 2. Expected recoil-electron spectra at SK for the best 
fitted parameters of the VEP scenarios, which are indicated 
in the legend of the plot as (sin^ 20g, I0A7I x 10^", 
preliminary data from SK are also plotted. 
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the best fitted parameters as well as the Xmin values for 
fixed and free fs- We found for fs = ^ that x^j^ = 1.49 
for 3-2=1 degree of freedom and for — 0.81 that 
Xmin = 0-32 for 3-3=0 degree of freedom. We also have 
checked that the allowed region for fixed *B flux (/s = 1) 
is rather similar to the one presented here and so we only 
give the values of the corresponding best fitted parame- 
ters for this case in Table |. 

Next we perform a spectral shape analysis fitting the 
^B spectrum measured by SK [|j using the following 
definition: 



E 



(4) 



where the sum is performed over all the 18 experimental 
points S°'"%Ei) normahzcd by BP98 SSM prediction for 
the recoil-electron energy Ei, ai is the total experimen- 
tal error and 5'"^''° is our theoretical prediction that was 
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TABLE I. The best fitted parameters and Xmm for the 
VEP induced long-wave length neutrino oscillation solution 
to the SNP. The local best fit points in the 2nd 90 % C.L. 
islands are indicated in the parentheses. 



Case 



sin^ 26g I0A7I X 10^ 



Rates (/s = 1) 1.0 (1.0) 1.71 (12.3) — 1.49 (1.88) 

Rates 1.0 (1.0) 1.70 (12.4) 0.81 (0.81) 0.32 (0.71) 

Spectrum 0.98 (0.39) 1.00 (0.24) 0.80 (0.66) 15.8 (19.8) 

Combined 1.0 (0.99) 1.65 (12.2) 0.82 (0.82) 22.0 (23.0) 

calculated using the BP98 SSM differential flux, the 
y — e scattering cross section , the survival probability 
as given by Eq. taking into account the eccentricity 
as we did for the rates, the experimental energy resolu- 
tion as in Ref. ||2^ and the detection efficiency as a step 
function with threshold i?th = 5.5 MeV. 

After the minimization with fs — 0.80 we have 
obtained Xmin ~ 1^.8 for 18-3 =15 degrees of freedom. 
The best fitted parameters that can be found in Table || 
permit us to compute the allowed region displayed in Fig. 
1 (b). 

Finally we have performed a combined fit of the rates 
and the spectrum obtaining the allowed region presented 
in Fig. 1 (c). Again we can read from Table | the best 
fitted parameters. We observe that the combined allowed 
region is essentially the same as the one obtained by the 
rates alone. In all cases presented in Figs. 1 (a)-(c) we 
have two isolated islands of 90% C. L. allowed regions. 
See Table | for the fitted values corresponding to the 
local minimum in these islands. We note that only the 
upper corner of the Fig. 1 (c), for |(/)A7| > 2 x 10^^"^ and 
maximal mixing in the v,, — > channel can be excluded 
by CCFR and moreover, there are no restrictions 
in the range of parameters we considered in the case of 
j/g — > i/^ or Vf, — > Us oscillations. 

In Fig. Ijwe show the expected recoil-electron spectrum 
in SK for various fitted parameters of the VEP solution 
to the SNP. We see that the data from the spectrum 
alone can be quite well described by the VEP oscillation 
mechanism (thick solid line), whereas the prediction for 
the best fitted parameters from the rates alone and from 
the combined fit give flatter curves (dashed and long- 
dashed lines). Nevertheless parameters for a "test point" 
taken inside the 90 % C. L. region of Fig. 1 (c) can give 
rise to some spectral distortion (thin solid line). 

We have performed the same analyses with rates as 
well as spectrum also for the — > Vs channel. Since 
the allowed regions as well as the fitted recoil-electron 
spectra obtained in this case are rather similar to the 
ones for active to active conversion, we do not present 
them here but only show the best fitted parameters and 
Xmin values in Table ||. Although the spectrum alone 
gives a comparable fit to the active to active case, we see 
that the rates can not be so well explained by this type of 
scenario and consequently the combination gives a worse 
fit. In spite of that this is still much better than the mass 
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FIG. 3. Yearly averaged survival probability for the best 
fitted parameters, indicated in the parentheses as (sin^2&(3, 
|(;/)A7| X 10^*) in the plots, which can explain well (a) the 
rates or (b) the SK spectrum by the VEP induced neutrino 
oscillation. The energies of the pp as well as ^Be neutrinos 
are also indicated by the thick dot-dashed and dashed line, 
respectively. 

induced active to sterile vacuum oscillation solution to 
SNP. 

To understand why it is possible to fit the solar neu- 
trino data we show in Fig. ^ (a) the survival probabilities 
for the best fitted parameters of the VEP induced oscilla- 
tion. Due to the specific energy dependence of the proba- 
bility assumed here we can actually strongly suppress the 
Be line and still keep the neutrino flux high enough 
to be in agreement with Ga data, and at the same time 
obtain 50 % reduction of the ^B neutrino flux, which 
is in fact the required suppression pattern of the solar 
neutrino fluxes in order to get a good fit [ p8[ . 

Because of the contributions from the strong smearing 
in energy of the scattered electron and of the finite exper- 
imental energy resolution, the probability alone can not 
give us a precise insight on the spectral shape. We can 
only qualitatively expect some distortion for the proba- 
bihty in Fig. | (b). 

TABLE II. Same as Table I but for the case of —> Va 
conversion. 



Case 



sin^ 2610 I0A7I X 10^ 



fB 



yimin 



Rates (/s = 1) 1.0 (1.0) 
Rates 1.0 (1.0) 

Spectrum 0.88 (0.33) 
Combined 1.0 (1.0) 



1.80 (12.1) — 3.06 (3.87) 

1.80 (12.1) 0.94 (0.94) 2.96 (3.85) 

1.01 (0.24) 0.84 (0.66) 15.6 (19.7) 

1.66 (12.5) 0.94 (0.94) 24.7 (26.2) 



Finally, let us discuss about the seasonal variation of 
the solar neutrino signal. In contrast to the usual vac- 
uum oscillation solution to the SNP, in this scenario, no 
strong seasonal effect is expected in any of the present or 
future experiments, even the ones that will be sensitive to 
^Be neutrinos such as Borexino and Hellaz j3^. Con- 
trary to the usual vacuum oscillation case, the oscillation 
length for the low energy pp and ^Be neutrinos are very 
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FIG. 4. Expected seasonal variations for the fitted pa- 
rameters of VEP scenarios, indicated in the parentheses as 
(sin^ 26'g, I0A7I X 10'^'', /s) in each plot. The preliminary 
data from SK are also plotted. Variations due to the eccen- 
tricity of the Earth orbit (~ were subtracted. 

large, comparable to or only a few times smaller than the 
Sun-Earth distance, so that the effect of the eccentricity 
in the oscillation probability is small. On the other hand, 
for higher energy neutrinos relevant for SK, the effect of 
the eccentricity in the probability could be large, but av- 
eraged out after the integration over a certain neutrino 
energy range. These observations are confirmed in Fig. 
^ where we present the expected seasonal variations for 
the best fitted parameters of the VEP induced oscillation 
scenario. 

In conclusion we found a new solution to the SNP 
which is comparable in quality of the fit to the other 
suggested ones. 
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